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p ,. ABSTRACT 

' ' We report the serendipitous discovery of a ring nebula around a candidate Wolf- 

\ Rayet (WR) star, HBHA 4202-22, in Cygnus using the Spitzer Space Telescope archival 

^— > i data. Our spectroscopic follow-up observations confirmed the WR nature of this star 

I (we named it WR138a) and showed that it belongs to the WN8-9h subtype. We 

thereby add a new example to the known sample of late WN stars with circumstellar 
nebulae. We analyzed the spectrum of WR 138a by using the Potsdam Wolf-Rayet 
(PoWR) model atmospheres, obtaining a stellar temperature of 40 kK. The stellar wind 
composition is dominated by helium with 20 per cent of hydrogen. The stellar spectrum 
is highly reddened and absorbed {Eb-v = 2.4 mag. Ay = 7.4 mag). Adopting a stellar 
I luminosity of logL/L© = 5.3, the star has a mass-loss rate of 10^"*'^ MQyr~^, and 

i resides in a distance of 4.2 kpc. We measured the proper motion for WR138a and 

found that it is a runaway star with a peculiar velocity of ~ 50 kms^^. Implications 
of the runaway nature of WR 138a for constraining the mass of its progenitor star and 
understanding the origin of its ring nebula are discussed. 
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1 INTRODUCTION method for searching for WR stars using optical surveys. In 

this respect, the infrared (IR) observations are of high im- 
It IS beueved that single stars or solar metallicity with an . ■ • j_ j_ n i- /toain • i 

° portance since the interstellar medium (IbM) is much more 

initial mass of > 20 — 25 Mo end their lives as Wolf-Rayet , ^ • m i -^i j_i i- i 

/, r , T o ir n trausparcut m the IR comparcd With the optical. 

(WR) stars (Vanbeveren, De Loore & Van Rensbergen 1998; 

Meynet & Maeder 2003). The relatively small number of The effectiveness of IR observations in identifying new 

massive stars and the short duration of the WR phase (< 0.5 ^R stars was demonstrated by Shara et al. (2009) who car- 

Myr) makes the WR stars rare objects. The present total "ed out a near-IR narrow-band imaging survey of the inner 

number of WR stars in the Galaxy does not exceed sev- Galactic plane to select candidate WR stars and discovered 

eral thousands (e.g. Shara et al. 1999; van der Hucht 2001). ^1 new WR stars in spectroscopic follow-ups (cf. Homeier 

The known population of WR stars is, however, an order of et al. 2003; Hadfield et al. 2007; Mauerhan, Van Dyk & 

magnitude less numerous (van der Hucht 2006). The large Morris 2009). Another possible way to search for new WR 

disparity between the expected and the observed number of (^n^ ot^i^r evolved massive) stars is through detection of 

WR stars is mostly caused by the huge interstellar extinction ^heir (IR) circumstellar nebulae, created after a massive star 

in the Galactic plane, which strongly Umits the traditional leaves the main sequence (e.g. Marston et al. 1999; Egan et 

al. 2002; Clark et al. 2003). 

In this paper, we report the discovery of a new WN8- 
* E-mail: vgvaram@mx.iki.rssi.ru (VVG); fabrika, olga, aza- 9^ star in Cygnus via the detection of its IR nebula. We 

mat@sao.ru (SF, OS, AFV); wrh, lida@astro.physik.uni- thereby add a new example to the known sample of late 

potsdam.de (WRH, LMO); goray, cher@sai.msu.ru (VPG, AMC); WN (WNL) stars with circumstellar nebulae (e.g. Esteban 

bomans@astro.ruhr-uni-bochum.dc (DJB) et al. 1993). Inspired by this discovery, we undertook an ex- 



2 V. V. Gvaramadze et al. 











_.ja: J 










■ 














■ 






























""C" — 














T 


































20.0 


10. 


3 


20:17:00.0 



Right ascension 

Figure 1. Spitzer MIPS 24 /^m image of a new ring nebula in 
Cygnus and its central star HBHA 4202-22 (WR138a). The ori- 
entation of the spectrograph slit in our optical observation is in- 
dicated by a line. 



tensive search for similar objects using the archival data of 
the Spitzer Legacy Program^J and discovered many dozens 
of ring-like and bipolar IR shells (Gvaramadze, Kniazev & 
Fabrika 2009; cf. Carey et al. 2009). Spectroscopic follow- 
ups of central stars associated with two dozens of the shells 
showed that they are either candidate Luminous Blue Vari- 
ables or related (WNL, blue supergiant) stars (Gvaramadze 
et al. 2009, in preparation; Kniazev et al. 2009, in prepara- 
tion) , which confirmed that the IR imaging provides a pow- 
erful tool for revealing evolved massive stars via detection 
of their circumstellar nebulae. 



2 A NEW RING NEBULA AND ITS CENTRAL 
STAR 

The nebula, which is the subject of this paper, was dis- 
covered serendipitously during our search for bow shocks 
around OB stars running away from star clusters in the 
Cygnus X region (for motivation of this search see Gvara- 
madze & Bomans 2008a). The nebula was detected in the 
archival data of the Spitzer Space Telescope, obtained with 
the Infrared Array Camera (IRAC) and the Multiband 
Imaging Photometer for Spitzer (MIPS) (Werner et al. 2004) 
within the framework of the Cygnus-X Spitzer Legacy Sur- 
veja. 

Fig.[T] shows the MIPS 24 /im image of the nebula and 
its central star. The nebula has a clear ring-like structure 
with a diameter of ~ 2.3 arcmin and enhanced brightness 
along the west rim. The central star is offset for ~ 0.2 ar- 
cmin from the geometric centre of the nebula, being closer 
to the brightest portion of the nebula (a possible origin of 
the nebula and its brightness asymmetry are discussed in 
Section[S3J . The Cygnus-X Spitzer Legacy Survey also pro- 



^ http:/ /irsa.ipac.caltech.cdu/Missions/spitzer.html 
^ http:/ /www. cfa.harvard.edu/cygnusX/indox. html 



vides IRAC 3.6, 5.8 ^m and MIPS 70 fim images of the field 
containing the nebula; none of them shows its signatures. 

The optical counterpart to the central star was identi- 
fied by Dolidze (1971) as a possible WR star due to the pres- 
ence in its spectrum of an emission band around A6750 A. 
This star was included in the Catalogue of H-alpha emission 
stars in the northern Milky Way (Kohoutek & Wehmeyer 
1999) as a candidate WR star (named in the SIMBAD 
database as [KW97] 44-47 or HBHA 4202-22), but was not 
listed in the Vllth Catalogue of Galactic Wolf-Rayet Stars 
(van der Hucht 2001) nor in the annex to this catalogue 
(van der Hucht 2006). Recently, Corradi et al. (2008) put 
this star in the list of candidate symbiotic stars on the basis 
of its IPHAS and 2MASS colours. 



3 OBSERVATIONS 

3.1 Spectroscopy 

To verify the nature of HBHA 4202-22, we had performed 
spectroscopic observations with the Russian 6-m telescope 
using the SCORPIc|f| focal reducer (Afanasiev & Moiseev 
2005) in a long-slit mode with a slit width of 1 arcsec, pro- 
viding a spectral resolution of 5 A. The orientation of the 
(6 arcmin long) slit is shown in Fig.[T] The spectra were 
taken on two occasions: 2008 October 10 (spectral range 
AA4000 - 5700 A) and 2008 November 11 (spectral ranges 
AA4000 - 5700 A and AA5700 - 7500 A), and reduced using 
standard procedures. 

3.2 Photometry 

Using secondary photometric standards in the regions of 
symbiotic variable stars V407 Cyg and AS 323 published by 
Henden & Honeycutt (1997), we determined accurate mag- 
nitudes of four comparison stars near HBHA 4202-22. The 
observations were carried out with the Special Astrophysical 
Observatory (SAO) 1-m Zeiss refiector and UBVRcIc pho- 
tometer in the good photometric night on 2009 May 29. The 
photometry accuracy is better than 0.02 mag. We have made 
expanded CCD standard star set of 18 stars to measure 
a DSS POSS-1 blue image. HBHA 4202-22 was measured 
using POSS-O plate taken in the epoch of 1954.57, CCD 
frames taken with the Russian 6-m telescope in 2008.90 and 
CCD frames taken with SAO 1-m Zeiss telescope in 2009. 
The results are given in Table[T] 

The uniform 1-m Zeiss photometry shows that the ob- 
ject is variable with amplitudes between 0.04 and 0.07 in dif- 
ferent filters in the time scale of a month. A total variability 
range in B filter between 1955 and 2009 is 0.19 mag. This 
level of variability is typical of WNL stars (e.g. Marchenko 
et al. 1998). 

To measure the Spizer's IRAC and MIPS fluxes from 
HBHA 4202-22, we performed aperture photometry of the 
star using the MOPEX/APEX source extraction package and 
applied aperture correction, obtained from the nearby bright 
point sources. The fluxes are listed in Table[2] the estimated 

^ Spectral Camera with Optical Reducer for Pho- 
tometrical and Interferometrical Observations; 

[http: / / www.sao.ru/hq/lsfvo/devices/scorpio/scorpiolitmll 
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Table 1. Photometry of HBHA 4202-22 (WRlSSa) 



JD 


B 


V 


Rc 


Source 


2434951 


17.15 ±0.06 






POSS-1 


2454796 


17.25 ±0.01 


15.44 ±0.01 




6-m/SCORPIO 


2454981 


17.34 ±0.02 


15.57 ±0.02 


13.90 ±0.02 


1-m Zeiss 


2455008 


17.27 ±0.02 


15.53 ±0.02 


13.85 ± 0.02 


1-m Zeiss 



Table 2. Details of HBHA 4202-22 (WRlSSa) 



RA(J2000) 


201^ 17™ 08? 12 


Dec.(J2000) 


41°07'27'.'0 


l,b 


78.3203, 3.1536 


B (mag) 


17.25 ± 0.01 


V (mag) 


15.44 ±0.01 


J (mag) 


10.15 ± 0.02 


H (mag) 


9.27 ±0.02 


Ks (mag) 


8.65 ±0.02 


[3.6] (mjy) 


168.0 ±5.0 


[5.8] (mJy) 


154.0 ±4.6 


[24] (mJy) 


45.7 ± 1.4 


[70] (mJy) 


< 28.1 



errors are ~ 2 — 3 per cent. For the flux at 70 ^m, we give 
a 3(T upper limit since HBHA 4202-22 was not detected at 
this wavelength. 

The details of HBHA 4202-22 are summarized in Ta- 
ble[2l The coordinates and the J, H, Ks magnitudes are 
taken from the 2MASS All-Sky Catalog of Point Sources 
(Skrutskie et al. 2006). 



3.3 Proper motion 

We measured the proper motion for HBHA 4202-22 using 
four epochs in the POSS-1 and POSS-2 Digitized Sky Sur- 
veys and the epoch of our 6-m/SCORPIO frames. The full 
epoch range is 54.3 years. 21 reference stars located within 
an angular radius of 2.5 arcmin from the object position 
were used. Typical rms deviations of reference stars in these 
frames are between 0.21 and 0.56 arcsec. Linear trends are 
well seen in both coordinates. To evaluate the stellar proper 
motion and its errors, we got a linear approximation of these 
trends using the method of least squires and derived mean 
rms deviations of stellar coordinates from these linear ap- 
proximations. The proper motion errors were estimated by 
dividing these mean deviations by the epoch difference. The 
result is Ha = — 4.9±1.2 masyr"^, — — 4.7±1.2 masyr"^. 

To convert the observed proper motion into the true 
tangential velocity of the star, we use the Galactic constants 
7?o = 8 kpc and 6 = 200 kms"^ (e.g. Reid 1993; Avedis- 
ova 2005) and the solar peculiar motion {Uq,Vq,Wq) = 
(10.00, 5.25, 7.17) kms-i (Dehnen & Binney 1998), and 
adopt a distance d = 4.2 kpc (see Section f44)) . We found 
that the star is moving in the west-southwest direction with 
a peculiar (transverse) velocity of ~ 50 ± 24kms~^ (here 
we assume that the errors are mainly due to the errors in 
our proper motion measurements). Note that the vector of 
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Figure 2. Blue (top) and red (bottom) spectra of HBHA 4202- 
22 (WR138a) with some lines indicated (upper tow). In the blue 
spectrum, the bottom row of vertical bars marks numerous He II 
lines while the middle one marks He I lines. In the red spectrum, 
the vertical bars mark He II lines, whereas very strong He I lines 
are not marked. Hi emissions are not marked on both spectra 
as well. CCD fringes are notable at wavelengths longer than 
6900 A. 

stellar peculiar velocity is oriented by chance almost exactly 
along the spectrograph slit (see Fig.[T]). 



4 A NEW WOLF-RAYET STAR - WR 138A 
4.1 Spectral type 

In Fig. [2] we present the spectrum of HBHA 4202-22 in the 
blue and red regions. The spectrum shows strong emission 
lines of hydrogen, He I and Hell. We detect emissions of 
Nil, Nlll, Cm, Niv, Si IV. No forbidden lines can be seen 
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in the spectrum. Many of the weaker hnes in the spectrum 
show P Cygni profiles, while the strongest lines are entirely 
in emission. The strong absorptions visible in the red spec- 
trum are telluric. Numerous diffuse interstellar bands (DlBs) 
are observed, in the blue region they are AA4428, 4726, 4762, 
in the red the strongest one is at 6280 A. 

Our spectra also show nebular emission lines along the 
whole length of the slit. Radial velocity in Ha line is con- 
stant along the slit and equal to — 2f ± 5 kms~^. From the 
Balmer decrement, we estimated the interstellar extinction 
and obtained Ay = 1.9 ± 0.1 mag (also constant along the 
slit), which is much less than Ay = 7.4 mag derived from the 
stellar spectrum (see Section 4.4). Hence, we conclude that 
the nebular line emission originates from the foreground and 
is not related to HBHA 4202-22 and its ring nebula. 

The dominance of helium and nitrogen emission lines 
indicates that HBHA 4202-22 belongs to the nitrogen (WN) 
sequence of WR stars. The presence of hydrogen emission 
lines is typical for WNL subtypes (Hamann et al. 1991). We 
named this star WR 138a, in accordance with the number- 
ing system of the VHth Catalogue of Galactic Wolf-Rayet 
Stars by van der Hucht (2001). To determine the subtype 
of WR 138a more precisely, we use the classification scheme 
proposed by Smith (1968) and updated for WNL stars by 
Crowther, HiUier & Smith (1995) and Smith, Shara & Moffat 
(1996) . The relative strengths of N ill A 4640, N iv A 4058 and 
He II A 4686 lines and the P Cygni profiles of He I lines sug- 
gest that the star belongs to the WN8-9 subtype. The WN8- 
9 classification for WR 138a also follows from the position of 
this star on the diagrams showing a comparison of the emis- 
sion equivalent widths (EW) of He I A5876 versus He ll A4686 
lines and Hell A4686 EW versus FWHM for WNL stars in 
the Milky Way and the Large Magellanic Cloud (LMC) (see 
Fig. 1 of Crowther & Smith 1997). With EW(5876)=40.0 A, 
EW(4686)=19.5 A and FWHM(4686)=7.7 A, WR138a lies 
almost exactly on the line separating the WN8 and WN9 
stars. The further evidence for this classification is given 
in Section H31 Following the three-dimensional classification 
for WN stars by Smith et al. (1996), we add a 'h' suffix to 
WN8-9 to indicate the presence of hydrogen emission lines, 
so that WR138a is a WN8-9h star. 

4.2 Radial velocity 

We have measured the radial velocities of the main emis- 
sion lines in the spectrum and the line widths corrected for 
the spectral resolution. We used the Gaussian analysis for 
the measurements. The accuracy of radial velocity measure- 
ments is < 10 kms~^, while the relative velocities are mea- 
sured with a notably better precision. 

For hydrogen lines H7, H/3 and Hq, we measured he- 
liocentric radial velocities Vr — —10 kms~^ (FWHM ~ 
610 ±30 kms"^). The strongest Hell A4686 line shows 1;^ ~ 
-37 kms"^ (FWHM ~ 490 kms"i) and it does not have 
the P Cyg profile. All He 11 lines with the P Cyg profile have 
Vr ~ 49 kms"^ (FWHM ~ 170 kms"^) and the radial ve- 
locity measured for the absorption peak Wa — —380 kms~^. 
Blue shifted absorption obviously shifts emission to the red 
side and makes it narrower. He I lines show the same be- 
haviour as Hen lines. The mean radial velocity for all 'non- 
PCyg' Hei lines is ~ -11 kms'^ (FWHM ~ 590 kms"^). 

We may adopt for the star's radial velocity Vr ~ 



—20 kms^ , which is an average over the hydrogen lines. 
Hell A4686 and the 'non-PCyg' He I lines. It is necessary 
to note, however, that the radial velocity of the star cannot 
be determined with good accuracy, because all the spectral 
lines axe formed in the wind. Resulting from a complicate 
formation process in the expanding atmosphere, the emer- 
gent line profiles are not exactly symmetric. Due to the high 
wind speed (700 kms~^, as found in our spectral analysis, 
see below), these effects are large compared to possible ra- 
dial velocity uncertainties. 

4.3 Spectral analysis and stellar parameters 

To analyze the stellar spectrum and to derive the funda- 
mental parameters of WR 138a, we use the Potsdam Wolf- 
Rayet (PoWR) models for expanding stellar atmospheres. 
These models account for complex model atoms including 
iron-line blanketing in non-LTE (for a detailed description 
see Hamann & Grafener 2004). For abundances of trace el- 
ements, we adopt values which are typical for Galactic WN 
stars - N: 0.015, C: 0.0001, Fe: 0.0014 (Hamann & Grafener 
2004). 

The main parameters of a WR-type atmospheres are 
the stellar temperature, T*, and the so-called transformed 
radius, Rt. The stellar temperature T* denotes the effective 
temperature related to the radius R*, i.e. L = aT^AnRl, 
where a is the Stefan-Boltzmann constant and 7?« is by def- 
inition at a Rosseland optical depth of 20. Rt is related to 
the mass-loss rate M and defined by 



* * 2500 kms-i / 10-4 Moyr-iJ 

where D is the clumping contrast (adopted here to be 4 
throughout the wind as a typical value for WN stars; see 
Hamann & Koesterke 1998), and Voo is the terminal velocity 
of the wind. 

These basic stellar parameters T, and Rt are deter- 
mined from fitting the lines in the normalized spectrum (see 
Fig-El • The normalization is in fact achieved by dividing the 
absolutely calibrated observed spectrum by the theoretical 
continuum, which makes the total procedure described in 
the following an iterative process. 

A first orientation about the proper choice of these pa- 
rameters can be obtained by comparing the observed spec- 
trum to the published grids of WN models (Hamann & 
Grafener 2004). For a detailed fit, we calculate individual 
models for this star with the adequate terminal wind veloc- 
ity, which is Voo = 700 kms~^ as inferred from fitting the 
width of the emission line profiles. 

The stellar temperature is adjusted such that the bal- 
ance between the lines from He I versus He 11 is reproduced. 
At the same time, Rt is adjusted, which infiuences the 
strength of the emission lines in general. The parameters of 
the best-fitting model are compiled in Table|3] (where for the 
completeness we also give the hydrogen ionising luminosity, 
c&i). Note that according to its location in the T* — Rt plane 
(Hamann & Grafener 2004), WR138a belongs to a spectral 
subtype later than WN8 (cf. Section|4TT} ■ 

Fig.|3]shows that the Balmer lines of hydrogen dominate 
their blend with lines of the Hell Pickering series. Radial 
velocities of the Balmer lines confirm this. These blends are 
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Figure 3. Observed optical spectrum (blue/solid line) of WR 138a, compared with the best fitting model (red/dotted line) with the 
parameters as given in Table[3] The absolutely calibrated observation had been divided by the reddened model continuum for normal- 
ization. 



nicely matched by models with a hydrogen mass fraction of 
20 per cent, leaving about 80 per cent of the mass for helium 
which is a typical composition for WNL stars. 

The spectral analysis alone cannot tell the absolute di- 
mensions of the star. These are related to the distance of the 
object, and will be discussed in the following subsection. 

4.4 Luminosity and distance of WR 138a 

The lowest plausible luminosity for a massive WNL star 
is about log L/Lq = 5.3 (Hamann, Grafener & Liermann 
2006). Let us adopt this value for a moment. This choice 
implies the values for the stellar radius and mass-loss rate as 
given in Table|31 Note that, in fair approximation, WR mod- 
els can be scaled to different luminosities, with _R» oc L^^^, 
M oc I/^^**, and d oc L^^^ (Schmutz, Hamann & Wessolowski 
1989; Hamann & Koesterke 1998). 

To estimate the corresponding distance to WR 138a, we 
now fit its observed spectral energy distribution (SED) with 



the model SED (Fig.|4ll. Two parameters can be adjusted 
for the fit, the distance d and the reddening Eb-v- The 
reddening law we adopt from Seaton (1979) in the optical 
and from Moneti et al. (2001) in the IR. A perfect fit to 
the whole SED is achieved with Eb-v ~ 2.4 mag, implying 
an extinction in the visual of Ay ~ 7.4 mag and d = 4.2 
kpc. The strong interstellar absorption also shows up in the 
observed spectrum by pronounced DIBs (see Fig. [2]). 

The line-of-sight towards WR 138a is nearly tangential 
to the local (Orion) spiral arm, whose extent in this direction 
is ~ 4-6 kpc (McCutcheon & Shuter 1970; Russeil 2003). A 
distance of about 4 kpc would imply that this star is located 
at the far end of the Orion arm, which is compatible with 
the high reddening. 

Alternatively, one might consider if our object is in fact 
not a massive star with ring nebula, but a planetary nebula 
(FN). FN central stars with WN-type spectra are very rare; 
only two of them are known in our Galaxy [FMR5, Morgan, 
Farker & Cohen (2003) and FB 8, Todt et al. (2009)] and 
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Table 3. Stellar parameters for WR 138a 



Spectral type 


WN8-9h 


T, [kK] 


40 


logilt [Rq] 


0.9 


Voo [kms^^] 


700 


log L[Lq] 


5.3 


R* [Re] 


9.4 


logM[M0yr-i] 


-4.7 


Eb-v [mag] 


2.4 


d [kpc] 


4.2 


log [s^-"-] 


48.9 




3.5 4.0 4.5 5.0 5.5 

log X [A] 



Figure 4. Observed flux distribution of WR 138a (blue/noisy) in 
absolute units, including the calibrated spectrum and the photo- 
metric measurements compiled in Table[2l compared to the emer- 
gent flux of the model continuum (red/smooth line), in the op- 
tical also shown with lines. The model flux has been reddened 
and scaled to the distance according to the parameters given in 
TableE] 



an eruptive one in the LMC [N66; Hamann et al. (2003)]. 
With a typical luminosity of a PN central star, 6000 Lq , the 
photometric distance would become only 0.5 kpc, placing the 
object into the foreground of the Orion arm. In this case, the 
strong reddening would have to be of circumstellar origin. 
Although we cannot strictly rule out this scenario, it appears 
very artificial and unlikely. 

The distance of d = 4.2 kpc, which relies on the adopted 
stellar luminosity of log L/Lq = 5.3, places WR138a at a 
height of 230 pc above the Galactic plane. As we will discuss 
in the next section, this is just compatible with the inference 
that WR 138a is a runaway star that has been ejected from 
the Galactic plane. Many WNL stars have a much larger 
luminosity than log L/Lq — 5.3 (Hamann et al. 2006). A 
higher luminosity, however, would increase the implied dis- 
tance (d oc L^^^) and hence lead to an even larger height 
above the Galactic plane, while on the other hand the evo- 
lutionary lifetime decreases with higher luminosity. Both ef- 
fects together tend to make the runaway scenario impossible 
for much larger d and L (see Sect. l5.l| l. Hence, we think that 
our choice of log L/L0 = 5.3 is the most plausible one. 



5 DISCUSSION 

5.1 WR 138a as a runaway star 

Our proper motion measurements for WR 138a (Section l3.3|l 
showed that the star is moving in the west-southwest di- 
rection (i.e. away from the Galactic plane) with a peculiar 
(transverse) velocity of ~ 50 kms~^ (for d = 4.2 kpc), which 
is typical of runaway stars. In the Galactic coordinates, the 
components of the peculiar velocity are almost equal to each 
other: vi ~ —36 kms^^, vi, ~ 34 kms~^. 

Assuming that the progenitor star of WR 138a was born 
near the Galactic plane and ejected from the parent star 
cluster soon after the birth, one has that the time it needs 
to reach the present height above the Galactic plane, ~ 
dsinb/vi,, should be ^ tn, where tn is the H-burning lifetime 
of the star. Adopting the minimum plausible luminosity for 
WNL stars of log{L / Lq) — 5.3, one finds the current mass 
of WR138a of ~ 13 Mq, which in turn implies the initial 
mass of the star of ~ 30 Mq and tn — 6.0 — 7.0 Myr (e.g. 
Meynet & Maeder 2003). Thus, one has that t,{~ 6.7 Myr) 
< tn, so that the star had enough time to reach its present 
location. 

Placing WR 138a at, for example, d = 6 kpc would re- 
sult in a twofold increase of the stellar luminosity and in a 
increase of the present and initial mass of the star to, re- 
spectively, ~ 17 Mq and ~ 40 Mq. In this case, i* ~ 7.0 
Myr > tn — 4.5 — 5.5 Myr. Further increase of d makes the 
discrepancy between t, and tn more severe. The discrepancy 
would be even more sever if the star spent several Myr in the 
parent cluster before it was ejected (Gvaramadze & Bomans 
2008b; Gvaramadze, Gualandris & Portegies Zwart 2009). 
From this follows that the most likely initial mass of the 
progenitor star of WR 138a was ~ 25 — 30 Mq. One cannot, 
however, exclude a possibility that the progenitor star was a 
blue straggler formed via merging of two less massive stars 
in the course of close binary-binary encounter, resulting in 
ejection of the merger product from the parent cluster (e.g. 
Gvaramadze & Bomans 2008b) . In this case, the distance to 
and the mass of WR 138a could be larger. 

The inference that WR 138a is a runaway star could 
also be used to understand the origin of its ring nebula. 

5.2 Origin of the nebula around WR 138a 

A runaway wind-blowing (e.g. WR) star moving through 
the ISM creates a highly elongated bubble (Weaver et al. 
1977; Brighenti & D'Ercole 1994), whose shape in the up- 
stream direction is determined by a ram pressure balance 
between the stellar wind and the ISM (Baranov, Krasnobaev 

6 Kulikovskii 1971). The characteristic scale of the leading 
edge of the bubble is Rwn — (AfwR«wR/5.67rmHnof*)^''^, 
where Mwr and wwr are the stellar mass-loss rate and 
the wind velocity, no is the number density of the ambi- 
ent medium, is the stellar peculiar velocity and ttih is 
the mass of a hydrogen atom. For Mwr ~ 10"* '^ MQyr~^ 
and uwR ~ 700 kms~^ (see Table[3]), — 50 kms~^ and 
no ~ 0.1 cm"3 (typical of the ISM at 2 = 230 pc; Dickey & 
Lockman 1990), one finds Rwr — 20 pc. The small size and 
the nearly circular shape of the nebula around WR 138a (at 
d = 4.2 kpc, the linear radius of the nebula Vneh — 1-4 pc) 
imply that the stellar wind interacts with the dense ambi- 
ent medium comoving with the star (cf. Lozinskaya 1992). 
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This consideration suggests that the immediate precursor of 
WR 138a was a red supergiant (RSG) star (i.e. the initial 
mass of the progenitor star was ^ 40 M©; cf. Section lSTTj) 
and that the WR wind still propagates through the region 
occupied by the dense material shed during the RSG phase. 

According to stellar evolutionary models, a massive star 
with the initial mass in the range from 25 to 40 M© evolves 
through the sequence O -> RSG ^ WN (e.g. Meynet & 
Maeder 2003). During the RSG phase, the star loses a con- 
siderable fraction of its initial mass in the form of slow, dense 
wind. The subsequent fast WR wind sweeps up the RSG 
wind and creates a circumstellar shell. For a spherically sym- 
metric RSG wind, the shell driven by the WR wind expands 
with a constant velocity (e.g. Chevalier & Imamura 1983) 
iigh (Mwri'wr^^RSg/3Mrsg)^''^, where Mrsg and wrsg 
are, respectively, the mass- loss rate and the wind velocity 
during the RSG phase. Adopting Mrsg = 3 x lO"'^ Moyr"^ 
and wrsg = 10 kms~^ (the figures typical of RSGs), one has 
Ush — 100 kms~^, which in turn gives us the dynamical age 
of the nebula rncb/wsh — 1.4 x 10* yr. The latter estimate 
suggests that WR 138a only recently entered the WR phase. 

The runaway nature of WR 138a allows us to propose 
a natural explanation of the brightness asymmetry of the 
nebula and the offset of the star towards the brightest por- 
tion of the nebula. Namely, we suggest that both are caused 
by the effect of the ram pressure of the ISM, which af- 
fects the zone occupied by the RSG wind by making it 
denser and more compact ahead of the star. This sugges- 
tion could be supported by an estimate of the upstream 
(minimum) size of the zone occupied by the RSG wind, 
i?RSG =i (MRSG«RSG/5.67rmHno«*)^''^ (~ 1.6 pc for the pa- 
rameters adopted above) ~ rneb, which shows that although 
the WR wind still propagates through the zone occupied by 
the RSG wind, it is already close to the edge of this zone 
to fill the density enhancement caused by the stellar mo- 
tion. An additional support for our suggestion comes from 
the proper motion measurements for WR 138a, which show 
that the star is moving in the direction implied by the above 
consideration. 

Our interpretation of the nebula around WR 138a as a 
circumstellar one (i.e. created via the wind-wind interaction) 
is consistent with the observational fact that small-scale 
(circumstellar) nebulae are predominantly associated with 
WNL stars (e.g. Lozinskaya & Tutukov 1981; Esteban et al. 
1993; Barniske, Oskinova & Hamann 2008), i.e. with young 
WR stars, whose winds interact with the material lost dur- 
ing the preceding evolutionary phases rather than with the 
ambient ISM. The scarcity of circumstellar nebulae around 
WR stars suggests that they are visible on a time-scale much 
shorter than the duration of the WR phase. It is conceivable 
to associate this time-scale with the crossing-time of the re- 
gion occupied by the RSG wind, i.e. ~ r-Rsa/vsh ^ 2— 5 x 10'' 
yr, where trsq = vrsg^rsg ~ 2 - 5 pc and t^SG ^ 0.2 - 0.5 
Myr is the duration of the RSG phase. 



6 SUMMARY 

We have serendipitously discovered a ring nebula around 
a candidate WR star, HBHA 4202-22, in Cygnus using the 
MIPS 24 ^xa data from the Spitzer Space Telescope archive. 
Our spectroscopic follow-up confirmed the WR nature of 



this star and showed that it belongs to the WN8-9h subtype 
(we named the star WR 138a). We analyzed the spectrum of 
WR 138a by using the Potsdam Wolf-Rayet (PoWR) model 
atmospheres, obtaining a stellar temperature of 40 kK. The 
stellar wind composition is dominated by helium with 20 per 
cent of hydrogen. The stellar spectrum is highly reddened 
and absorbed {Eb-v ~ 2.4 mag. Ay = 7.4 mag). Adopt- 
ing a stellar luminosity of log L/L© = 5.3, the star has a 
mass-loss rate of 10"*'^ M0yr~'^, and resides in a distance 
of 4.2 kpc. We measured the proper motion for WR 138a 
and found that it is a runaway star with a peculiar velocity 
of ~ 50 kms~^. The runaway nature of WR138a was used 
to constrain the initial mass of its progenitor star and to 
propose an explanation of the origin of its circular nebula. 
We found that the most likely initial mass of the progen- 
itor star was 25 — 30 M© (i.e. the immediate precursor of 
WR 138a was a RSG star) and suggested that the nebula 
around WR138a is a circumstellar one, created via the in- 
teraction between the WR wind and the dense material shed 
during the preceding RSG phase. 
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